#include <iostream>
#include "interpolation_polynomial.h"
#include <vector>
using namespace std;

int main()
{
  vector<double> x = {-5,-20.0/6,-10.0/6,0,10.0/6,20.0/6,5};
  vector<double> y = {1.0/26,9.0/109,9.0/34,1,9.0/34,9.0/109,1.0/26};
  Newton_formula p(x,y);
  vector<vector<double>> f = p.getf();
  for(int i = 0; i < x.size();i++)
    {
      cout<< f[i][i]<< " ";
      cout<< endl;
    }
      cout<< endl;
  return 0;
}
